Subthreshold pair production in short laser pulses 
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Abstract 

The e + e~ pair production by a probe photon traversing a linearly polarized laser pulse is treated 
as generalized nonlinear Breit- Wheeler process. For short laser pulses with very few oscillations of 
the electromagnetic field we find below the perturbative weak-field threshold y/s = 2m a similar 
enhancement of the pair production rate as for circular polarization. The strong subthreshold 
enhancement is traced back to the finite bandwidth of the laser pulse. A folding model is developed 
which accounts for the interplay of the frequency spectrum and the intensity distribution in the 
course of the pulse. 
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I. INTRODUCTION 



The Breit- Wheeler process [1], 7' + 7 — > e + + e~, describes pair production as a t channel 
process, which is sometimes termed "conversion of light into matter" . Considering it as 2-to-2 
process in a particle physics language, the Breit- Wheeler pair creation is a threshold process, 
i.e. the available energy, expressed by the Mandelstam variable s = (k' + k) 2 = 2u'u(l — 
cos0£,g), must fulfill s > 4m 2 , where m is the electron mass. Supposed, 7 describes optical 
(laser) photons with energy w ~ 1 eV, the energy of the counter propagating photon 7' must 
obey uj' > 250 GeV. While such high energy photons might be generated by various processes 
at high energy accelerators/colliders, e.g. in the time reversed Breit- Wheeler process (pair 
annihilation) or Compton backscattering of a laser beam off a very energetic electron beam, 
or in astrophysical environments, the availability of laboratory based 250 GeV photon sources 
is quite scarce. Nevertheless, the Breit-Wheeler process has been identified experimentally. 
The SLAC experiment E-144 [2] was a special set-up of one of the above mentioned options. 
In a sufficiently strong laser field, however, multi-photon processes are enabled [HHH], 
schematically 7' + wy — > e + + e~, often termed nonlinear Breit-Wheeler process. Indeed, 
the interpretation [7] of the SLAC experiment E-144 has proved this possibility. In fact, 
at least n = 4 laser photons were needed to produce a pair. For an all-optical set-up with 
oj' ~ u, n > 10 11 laser photons would be required when replacing in the above threshold 
formula k — > nk. Strong laser fields are necessary for such nonlinear effects. The electron 
and positron in a strong laser field acquire an effective mass, m* = m^Jl + Oq/2, where ao 
is the dimensionless laser-strength parameter (cf. [BI-fTO"] for a recent analysis of the mass 
dressing in strong pulsed laser fields). That is, the above threshold estimate would read 
s > 4m. 2 in disfavor of achieving the threshold for strong optical laser fields. The up-shift 
of the threshold, however, is compensated to some extent by the higher harmonics with 
thresholds s n = Am 2 /n. That means, in principle, the higher harmonics (labelled by n) 
seem to shift the threshold towards arbitrarily small values of s. For a given value of s, 
the minimum number of laser photons no is the smallest integer such that s no < s. In the 
perturbative regime, where ao <C 1, the harmonics are suppressed by factors of a^ n , and a 
noticeable pair production therefore extents not too far below the 2-to-2 threshold of 4m 2 . 
In the deep subthreshold region ao ^ 1 and k < 1, a large number of laser photons partic- 
ipate in the formation of pair with cross section a oc e -8 ^ 3 "^, where k = sag /(2m 2 ) is the 
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Ritus parameter for pair production [5]. It has a similar non-analytic dependence on the 
field strength parameter encoded in k as the Schwinger effect which depends on the electric 
field amplitude (cf. [EE].) At fixed ao, there is a steep decline of the pair creation probability 
for s < 2m 2 /ao, i.e. large values of ao shift the region of noticeable pair production to small 
values of s. 

The cross section for pair production in the nonlinear Breit- Wheeler process by an unpo- 
larized probe photon 7' in a plane wave with linear polarization reads in such a case [5] 

v(s)=J2n>n Vn(s) with 

a n (s) = % P&p f " {Al + a 2 (2u - 1) (a\ - A A 2 )\ , (1) 

where Aj = Aj(n, a,b) = ^- dt cos- 7 1 exp{i(a sin t — b sin 2t — nt)} are generalized Bessel 



functions with arguments n, a = ^^aoy 1 + a\j2 yu(u* n — u) cos<^, and b = a\umj 's; u* n is 
defined by u* n = n-^. Note that the relevant energy variable is still s = (k' + k) 2 . In Fig. 1, 
a survey of the cross section is exhibited as a contour plot over the ^/s vs. ao plane. The 
apparent structures visible in the region 0.8 < i/i/MeV < 1 for a < 1 are related to the 
turn over from the first to the second harmonic. The features propagate to increasing values 
of y/s for a > 1. The imprints of the onset of the second and third harmonics are still 
visible. In the weak-field region, ao < 1, the steep decline of the cross section as a function 
of \fs is evident below = 2m . 

In the case of circular polarization, Eq. ([T]) becomes modified [6]. The overall pattern 
remains, but with less pronounced structures due to azimuthal symmetry. In the strong-field 
region, a ^> 1, the subthreshold production with a large number of photons is suppressed 
due to the large angular momentum of the pair [5] . 

The hitherto discussed examples are for infinitely extended plane waves. In reality, strong 
laser fields are generated presently by the chirped pulse amplification technique, i.e. pulse 
compression, and the asymptotic final state refers to a free electron-positron pair, where m* 
does not matter. Therefore, a modified threshold behavior is to be expected. Moreover, 
a pulse of finite duration is not longer monochromatic, instead the power spectrum has a 
support of finite width. Also this effect will have imprints on the pair production, as does 
the variation of the intensity in the course of the pulse. Considering the (nonlinear) Breit- 
Wheeler process as cross channel of the (nonlinear) Compton scattering one can expect 
similar strong effects of the temporal pulse shape, as has been found for the latter one [12) - 




10" 1 



10° 



10 1 



FIG. 1: Contour plot of the cross section ([I]) over the y/s vs. ao plane. The numbers label curves 
of constant cross sections in mb. Dashed curves depict the thresholds s n for n = 1 • • • 5 (from top 
to bottom). 

ITS] . In fact, the differential e + e~ spectra in the (general) Breit- Wheeler process depend 
sensitively on the laser pulse shape j8]. Even for weak laser intensities ao < 1, a significant 
enhancement of the total pair production rate just below the threshold \fs = 2m has been 
found recently [16] for short circularly polarized pulses. 




Given this motivation we study here the pair production off a probe photon in a linearly 
polarized laser pulse. In section 2, we recapitulate the basic formulas within the Furry 
picture, where the interaction with the probe photon 7' is treated perturbatively, while the 
interaction of the electron and positron with the laser pulse is accounted for by Volkov states. 
The numerical results are discussed in section 3 for the total pair production probability. 
A folding model is introduced in section 4 to access to the effect of enhanced subthreshold 
pair production, thus identifying the relevant physical mechanisms. The summary is given 
in section 5. 
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II. PAIR PRODUCTION IN PULSED LASER FIELDS 



The pair creation process is a crossing channel of the Compton scattering. Accordingly, we 
evaluate the cross section of the process as decay of a probe photon 7' with four-momentum 
k' and polarization four- vector e' into a pair e + (7) + e~(7), where e ± ( , y) are laser dressed 
Volkov states (cf. p2]) which encode the interaction with the laser field. For weak laser 
fields, ao <C 1, this process may be resolved perturbatively into a series of diagrams, where, 
in addition to the probe photon 7', n laser photons 7 interact with the outgoing e . The 
diagram with n = 1 at a C 1 reproduces the Breit- Wheeler process. 
The linearly polarized laser pulse is described here by the real four-vector potential 

^(0) = ^e^(0)cos0 (2) 

with transverse polarization four- vector e M obeying k-e = 0; e denotes the elementary charge. 
The envelope function g{4>), with = k ■ x as invariant phase, encodes the temporal shape 
of the pulse with wave four-vector k and normalization g(0) = 1. We chose in this paper, 
following [8], 

gift = cos 2 MM , (3) 

for |0| < ttN and g(<p) = for |0| > ttN, where N is number of cycles in the pulse. In 
line with the notation in [16] we denote the considered process as finite pulse approximation 
(FPA), in contrast to the infinite pulse approximation (IPA) where g — > 1. Both, IPA and 
FPA, use plane waves, i.e. e • k — 0. 

The S matrix for the generalized nonlinear Breit- Wheeler process is 

d 4 xfy(x)^-^=^_ p (x), (4) 

where and denote the Volkov wave functions (to be taken with the potential of 
the outgoing electron and positron with four-momenta p' and p, respectively. For instance, 
the positron wave function is given by 

*_ p (x) = - 2^4(0)) (5) 

1 * 

= 2-7^/ #'[2ep • Atf) + e 2 A\<p')} (6) 


with outgoing free-field spinor v p . Feynman's slash notation is used, e.g. 4- = l^A 1 = 7 • A, 
where 7 stands for the Dirac matrices. We employ the light cone coordinates x± = x° ± x 3 



and 2 ) and orient the coordinate system to have k ~ (u, 0, 0, —u), i.e. k- = 2oj, 

k +; ± = and <p = k_x + /2. Therefore, the integration measure is d 4 x = dx_d 2 x±d<j)/k_. 
After integrating over the components x_ and xj_ the matrix element reads 



S = . lC (2tt) 4 / ^ 4 \p + p' -k' -£k)M(£), (7) 
^2k' 2p 2p' J 2tt 

M(£) = u p ,fv p B (£)+ma u p {^ - ^L) VpBl ({) ~ ajfeT^ vH^ffl (») 

with the three phase integrals 

B m (£) = J d<Pg m ^) cos m (0) exp {i (£0 - + . (9) 

The integral -Bo(^) does not contain an envelope function in the preexponential and needs 
a regularization prescription, e.g. as proposed in [12J. The energy-momentum conservation 
in ([7]) fixes the value £ = £q — {p~ + P- ~ k'_ ) / . The pair emission probability reads with 
a = e 2 /(4vr) 

dH/ = m^w^V) lM]2p ^ dvdv - (10) 

Here, we employ the rapidity y = |ln(p° +p 3 )/(p° — p 3 ) = |lnp + /p_ = In p + /m± and 
transverse momentum p± = Jp\ + p\ to parameterize the final state phase space of the 



positron; m± = ym 2 + p\ is the transverse mass. Due to the finite laser pulse in FPA, the 
variables p± and y are independent. That means, the three variables p±, y and I are related 
by m±(y,£) = £ e v^-y • m IPA (cf. [5]), i.e. for g — > 1, due to the periodicity of the laser 
pulse, the phase integrals B m collapse to a 5 comb and J ^ in (jT]) changes into a sum over 
harmonics n with £ — > i n = n — ^m 2 a^/ (k ■ p) — l,m 2 a^/ (k ■ p'). Inserting £ n into £ n ) 
and solving for p± one recovers the usual harmonics in IPA. 

III. TOTAL CROSS SECTION 



The numerical evaluation of the probability Eq. (10) leads to the cross section a = W/ (Jp 7 ), 
.1 = 2 and p 7 = m 2 a 2 ) /(27ia)A(f) = uj~ 2 J^dcftT 00 , where T 00 is the energy density of the 
laser background field, thus defining an equivalent effective laser pulse duration A(f>. For not 
too short pulses, where the approximations of [16] are suitable, one finds A0 = Jf^ d(pg 2 (4>). 
These definitions ensure that the cross section is always based on the same energy in the 
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FIG. 2: Contour plot of the cross section for linear polarization over the \fs vs. ao plane. The 
numbers label curves of constant cross sections in mb. FPA with N = 1. 

laser pulse irrespectively of its duration. The cross section is exhibited in Fig. [2] for iV = 1, 
i.e., only one cycle in the pulse, as contour plot over the \fs vs. ao plane. 
The white region, where a < 1 mb, is completely changed in FPA (Fig. [2]) in comparison 
with the IPA case (Fig. [T]). For a < 1, significant pair production occurs at smaller values 
of y/s. The reason is explained in [TH]. In IPA, the nth harmonic contributes only for s > s n ; 
in FPA, it has a finite contribution also at s < s n . This effect is most dramatic for small 
values of ao, where only the lowest harmonics contribute. 

To show the influence of the number of cycles in the pulse, we exhibit in Fig. [3] the cross 
sections u{\fs) for N = 1,3 and 5 and compare with the IPA case for two values of ao- 
The enhancement of the subthreshold production can be many orders of magnitude above 
the IPA value at small values of ao, in particular for very short pulses. The enhancement 
is stronger than that found in [16] for circular polarization, e.g. for ao = 0.01, our N = 5 
curve exceeds the second IPA harmonic by a factor of 20 at a/s = 0.9 MeV (see left panel 
of Fig. |3|, while in [T6] the A^ = 5 curve is even slightly below the IPA curve at the same 
value of y/s. (This can be attributed mainly to the different envelope function in [IS]-) For 
ao > 1, the pulse length does not matter for the chosen pulse shape: The differences of IPA 
and FPA are within a factor of two. 
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FIG. 3: Cross sections as a function of y/s for ao = 0.01 (left) and 0.1 (right). The dotted curves 
are for IPA, while the other curves are for FPA with N = 1 (solid red), 3 (dashed blue) and 5 
(dot-dashed green). 

IV. FOLDING THE POWER SPECTRUM WITH HARMONICS 

The above discussed enhancement of subthreshold pair production for a 1 is explained 
for a circularly polarized laser beam in [TH] by exploiting properties of the special functions 
entering the expressions for the S matrix element or cross section, respectively. Instead 
of transferring the formalism of [IE] to linear polarization, we apply here a simple folding 
model to explain in an alternative and more intuitive manner the effect. The key is the 
finite bandwidth of a pulse: The power spectrum contains frequencies lower and larger than 
the central value u. This is accounted for in our model by calculating the weighted average 
of the pair production cross section of all energy components which are present in the laser 
pulse. Thus, a laser pulse with larger bandwidth will lead to an increased subthreshold 
production due to the high-frequency content. Denoting by G(£) the Fourier transform of 
the pulse envelope, G(£) = J d<f)g(<j)) e ^, the weighted cross section for the nth harmonic 
can be defined as 

with the IPA cross sections in the weak-field (perturbative) limit, i.e., the leading term 
in an expansion in powers of ao- For instance, is the Breit- Wheeler cross section (cf. |17j ) 



(7 ( 0)(s) = (3 _ ^ ln _ 2A(2 _ ^ ^ (12) 
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FIG. 4: Comparison of the folding model Eq. (11) (solid curves for a = (a\) + {02), i.e. only 



the first two harmonics are taken into account) and the full numerical calculations from section 
III (dashed curves). The black dotted curves depict the modified nonlinear Breit- Wheeler cross 
sections R n for n = 1 and 2 (upper and lower curves). Left panel: ao = 0.01, right panel: 
a = 0.1. 



and the second harmonic reads for linear polarization [5] 
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where u n = ns/(Am 2 ) and (3 n = Jl—u~ l . In (11), R n = f g 2n ((fi)d(j)/ J g 2 ((p)d(f) is an 
average over the intensity of the pulse relevant for the nth harmonic. It depends on the 
pulse shape but is independent of the pulse length. (For instance, for the pulse ([3]), we find 
i?2 — 0.72916 and -R3 = 0.6016, while for a Gaussian pulse shape one has R n = the 
cosh pulse of [IB] is characterized by R 2 = 2/3 and R3 — 8/15; R\ — 1 by definition.) That 
means, the relative importance of different harmonics depends on the shape of the pulse. If 
the frequency distribution does not matter, such as for long pulses with N ^> I, the effect 



of the intensity distribution for ao C 1 causes (cr n )(s) — > R n ^^ in (11). The effect of 



the intensity distribution in the course of a pulse leads to an overall reduced subthreshold 
production, while the frequency spectrum can lead to an increased subthreshold production 
in the vicinity of the corresponding thresholds s n . This explains why FPA is below IPA for 
certain values of ^/s in Fig. [3] (e.g., at ^J~s ~ 0.8 MeV in the right panel for iV = 3 and 5). 
The comparison of that folding model with the above full numerical calculations is exhibited 
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FIG. 5 1 CTi ELS cl function of yfs for various values of A?" (solid curves). The dashed curve depicts 
(<7i)int for ao = 1. The dotted curves are the IPA results for the first harmonic ("IPAl": ao = 1, 
"IPAO": a 0). 

in Fig. |4| Despite of the simplicity of the folding model, the quantitative agreement is 
surprisingly good. We assign the deviations at small values of y/s to the omission of higher 
harmonics. 



In the model (11), we employed the independence of the intensity variation and energy 
distribution (frequency bandwidth) for small values of ao and also omitted effects arising 
from the shifted threshold 2m*. Although a pair emerges finally with asymptotically free 
e masses m, the mass shift during the production of the pair still leaves imprints on the 
differential spectra in a pulsed field [S]. Focussing now on these intensity effects for larger 
values of a , one can define (o"i) int = / d<fi g 2 (0)cri(s, a —> a g((p))/ J d<pg 2 (0) , where ai(s, a ) 
is the non-perturbative cross section in IPA related to the first harmonic, i.e. o\ from Eq. ([I]). 
That folding model approximates quite well the spectrum obtained in the previous section. 
For a comparison with numerical results, a suitable cut in phase space, e.g. i < 1.2 (this 
value is dictated by the lowest gap in the p± — y differential spectrum for not too large 
values of ao and not too short pulses), separates here the first FPA harmonic thus defining 
In fact, fig. [5] clearly exhibits the shift of the threshold from 2m* towards 2m even for 
ao = 1, i.e., the IPA threshold 2m* is replaced by the FPA threshold 2m for smooth pulses. 
However, already at the IPA threshold y/s = 2m* the cross section (fTi)mt starts dropping 
when going to smaller values of y/s. This behavior is much more prominent for pulses with 
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a pronounced flat-top section. The shape of the curve (01) 

int 

function of a/s again is 

sensitive to the pulse shape but not to N. The excess at yfs < 2m is due to effects of the 
frequency spectrum discussed above. 

This simple model consideration demonstrates the physics content of the pair production 
process in pulsed laser fields by identifying two regimes (i) N^ 1 > Oq, where the frequency 
distribution essentially determines the spectrum, and (ii) N^ 1 < al, > 1, where the 
intensity variation in the pulse plays an important role. Our findings may help to improve 
simulations of cascade formation [TS] which are based on certain parameterizations of the 
pair production cross section. 

V. SUMMARY 

In summary we have shown that for short weak and medium-intense laser pulses the gen- 
eralized nonlinear Breit- Wheeler process is strongly (order of magnitudes) enhanced in the 
subthreshold region a/s < 2m. The effect is similar for linear and circular polarizations. 
While a qualitative explanation was given in [16] for circular polarization, one may attribute 
the enhanced subthreshold pair production to the non-monochromaticity of the laser pulse 
due to its finite duration. This complements the strong impact of the temporal pulse shape 
on the multi-differential spectra found in [SJ. Effectively, the high-frequency content of a 
pulsed laser field leads to a decrease of the threshold harmonic no for given value of a/s. 
The nonlinear (multi-photon) Compton process as crossing channel of the Breit- Wheeler 
process has shown an analog sensitivity to the temporal beam shape [121415] . Our analysis 
evidences that pair production in pulsed laser fields depends sensitively on both the fre- 
quency spectrum and the intensity variation in the course of the pulse. Modelling the pulse 
by a box shape one would ignore the intensity variation which is particularly important for 
stronger laser fields with ao 3> 1 and higher harmonics with n> 1. The pair production 
is an important step in the formation of QED cascades [18] which are thought to become 
important in high-intense laser fields. 

The authors gratefully acknowledge stimulating discussions with T.Heinzl, A. Hosaka, 
H. Ruhl, and H. Takabe. We thank T. E. Cowan and R. Sauerbrey for their supportive 



11 



interest. 



[1] J. A. Wheeler, G. Breit, Phys. Rev. 46 (1934) 1087. 
[2] D. L. Burke et al., Phys. Rev. Lett. 79 (1997) 1626; 

C. Bamber et al., Phys. Rev. D 30 (1999) 092004. 
[3] H. R. Reiss, J. Math. Phys. 3 (1962) 59. 

[4] N. B. Narozhny, A. I. Nikishov, V. I. Ritus, Sov. Phys. JETP 20 (1965) 622; 
[5] V. I. Ritus, J. Sov. Laser Res. (United States), 6 (1985) 497. 
[6] D. Yu. Ivanov, G. L. Kotkin, V. G. Serbo, Eur. Phys. J. C 40 (2005) 27. 
[7] H. Hu, C. Miiller, C. H. Keitel, Phys. Rev. Lett. 105 (2010) 080401; 

A. Ilderton, Phys. Rev. Lett. 106 (2011) 020404. 
[8] T. Heinzl, A. Ilderton, M. Marklund, Phys. Lett. B 692 (2010) 250. 
[9] F. Hebenstreit, A. Ilderton, M. Marklund, J. Zamanian, Phys. Rev. D 83 (2011) 065007. 
[10] A. Ilderton, Few Body Syst. 52 (2012) 431; 



C. Harvey, T. Heinzl, A. Ilderton, M. Marklund, arXiv:1203.6077 [hep-ph]. 
[11] A. Ringwald, Phys. Lett. B 510 (2001) 107. 

[12] M. Boca, V. Florescu, Phys. Rev. A 80 (2009) 053403. 

[13] F. Mackenroth, A. Di Piazza, Phys. Rev. A 83 (2011) 032106. 

[14] N. B. Narozhnyi, M. S. Fofanov, J. Exp. Theor. Phys. 83 (1996) 14 [Zh. Eksp. Teor. Fiz. 110 
(1996) 26]. 

[15] T. Heinzl, D. Seipt, B. Kampfer, Phys. Rev. A 81 (2010) 022125; 

D. Seipt, B. Kampfer, Phys. Rev. A 83 (2011) 022101. 



[16] A. I. Titov, H. Takabe, B. Kampfer, A. Hosaka, arXiv: 1205.3880 [hep-ph], Phys. Rev. Lett., 
in print. 

[17] L. D. Landau, E. M. Lifshitz, Course of Theoretical Physics, vol. 4, Akademie Verlag, Berlin, 
1980. 

[18] A. R. Bell, J. G. Kirk, Phys. Rev. Lett. 101 (2008) 200403; 

A. M. Fedotov, N. B. Narozhny, G. Mourou, G. Korn, Phys. Rev. Lett. 105 (2010) 080402; 
E. N. Nerush, I. Y. Kostyukov, A. M. Fedotov, N. B. Narozhny N. V. Elkina, H. Ruhl, Phys. 
Rev. Lett. 106 (2011) 035001; 

12 



N. V. Elkina, A. M. Fedotov, I. Y. Kostyukov, M. V. Legkov, N. B. Narozhny, E. N. Nerush, 
H. Ruhl, Phys. Rev. ST AB 14 (2011) 054401. 



13 



